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Abstract 

The matrix elements for K — > 'kttIv decays are described by four form fac- 
tors F, G, H and R. We complete previous calculations by evaluating R at 
next-to-leading order in the low-energy expansion. We then estimate higher 
order contributions using dispersion relations and determine the low-energy 
constants L\, L2 and L3 from data on K e $ decays and on elastic pion scatter- 
ing. Finally, we present predictions for the slope of the form factor G and for 
total decay rates. 



(J Work supported in part by Schweizerischer Nationalfonds/ Bundesamt fur 
Bildung und Wissenschaft (BBW)/ EEC Human Capital and Mobility Pro- 
gram/ Fondazione "Angelo della Riccia"/ INFN. 

a) NORDITA, Blegdamsvej 17, DK-2100 Copenhagen, Denmark. 

b) Universitat Bern, Sidlerstrasse 5, CH— 3012 Bern, Switzerland. 

c) Dipartimento di Fisica, Universita di Roma II - "Tor Vergata", Via della 
Ricerca Scientifica 1, 1-00173 Roma, Italy. 



1 Introduction 



In this article we analyze K14 decays, 

K — > inxlv ; I = e, /1 , (1-1) 

in the framework of chiral perturbation theory (CHPT). This method, also called 
"energy expansion" in the following, is based on an expansion of the Green functions 
in powers of the external momenta and of the light quark masses The matrix 

elements for processes are described by four form factors F, G, H and R. Their 
energy expansion reads 

I= M_K | /( 0) +/ (2) +J (4) + ...| ;I = FjG}HjR , (1.2) 

where is a quantity of order E n . The predictions for the lowest order terms were 
first given by Weinberg ||. The anomaly contribution has been determined in 
0, whereas F^ and have been evaluated in || [J. For a calculation of 
see [[TTJ. The experimental results |TTJ for F, G turn out to be 30 — 50% above the 



leading contributions. The missing piece must then come from higher orders. The 
expressions for F^ 2 \ G^ involve the low-energy constants L±, . . . , Lg, which are not 
fixed by chiral symmetry alone and which must be determined phenomenologically 
at the present stage of our capability to solve QCD. This was done in ||, where 
the Lj's were pinned down using experimental data (not related to decays) 
and involving the large- Nq prediction which states that, in the limit where the 
number of colors becomes large, certain combinations of the low-energy constants 
are suppressed. The decay ( |1.1|) is the simplest process where this rule can be tested 
H, 0]. In addition, it allows one to perform an independent determination of Li, L2 
and L 3 and thus to check consistency with other data. 

The aim of the present article is threefold. First, we fill the gap in the literature 
and evaluate also the next-to- leading order term R( 2 \ [The amplitude R is com- 
pletely negligible in K ei decays, because its contribution to the rate is suppressed by 
the factor mf. It must be retained, however, in the channel ]. Second, we note 
that, because the strange quark mass is not very small on a typical hadronic scale, 
the corrections 1^ to the leading-order terms of the form factors are substantial. 
The determination of the Lj's from decays is therefore affected with substantial 
uncertainties if carried out using only the first two terms in the expansion ( p..2[ ), 
as was done in || || . Here we improve these calculations by estimating the size of 
higher order contributions to F . We use for this purpose the method developed in 
|T2H , which amounts to write a dispersive representation for the relevant partial wave 
amplitudes, fixing the corresponding subtraction constants with chiral perturbation 
theory. We are then able to reduce the uncertainties in the determination of L 1; L 2 
and L3 in a significant manner, even more so if data on elastic hit scattering is 
considered in addition. Third, we predict the slope of the form factor G, and show 
that we may evaluate total decay rates for all channels in decays within rather 
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small uncertainties, provided the leading S- and P- wave form factors have been 
determined experimentally e.g. from K + — > 7r + 7r _ e + i/ e decays. 

The plan of the paper is as follows. In section 2, we provide the necessary 
kinematics and the definition of the form factors. Section 3 contains the result of 
the one-loop calculation of these quantities. In section 4, we use dispersion relations 
to construct a I — 5*-wave amplitude which has the correct phase to higher 
orders in the low-energy expansion. In section 5, we use this improved amplitude to 
determine the low-energy constants Li,L 2 and L 3 . Section 6 contains predictions, 
whereas a summary and concluding remarks are presented in section 7. 

2 Kinematics and form factors 

We discuss the decays 

K+(jp) - n+( Pl ) 7r(p 2 ) l + ( Pl ) v^vv) , (2.1) 
K+(p) - ttV) *°(P2) 1 + (Pi) m(j?u) , (2.2) 
K°{p) -> 7T°(pi) 7r-(p 2 ) Z+(p,) utfa) ; Z = e,/x , (2.3) 

and their charge conjugate modes. We do not consider isospin violating contributions 
and correspondingly set m u = m d , a QED = 0. 

2.1 Kinematics 

We start with the process ( |2.1| ). The full kinematics of this decay requires five 
variables. We will use the ones introduced by Cabibbo and Maksymowicz [fL3| . It is 
convenient to consider three reference frames, namely the K + rest system (Er-), the 
7r + 7r~ center-of-mass system (S 27r ) and the l + vi center-of-mass system (S^). Then 
the variables are 

1. s n , the effective mass squared of the dipion system, 

2. si, the effective mass squared of the dilepton system, 

3. the angle of the ir + in E 27r with respect to the dipion line of flight in E^, 

4. 9i, the angle of the / + in Ej w with respect to the dilepton line of flight in E^, 
and 

5. 0, the angle between the plane formed by the pions in E^ and the correspond- 
ing plane formed by the dileptons. 

The angles 9 W , 9i and are displayed in Fig. 1. 
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The range of the variables is 

4M 2 < SlT = { P i+p 2 ) 2 < {M K - mi ) 2 , 
m] < si = {pi + p u f < (M K - , 
< 9„, 9i < 7T, < cj) < 2?r. (2.4) 

It is useful to furthermore introduce the following combinations of four vectors 

P = Pi+P2, Q = Pi-P2, L = pi+p„, N = pi-p u . (2.5) 

Below we will also use the variables 

t= (p!-p) 2 ,u = (p 2 -p) 2 ,v = t-u. (2.6) 

These are related to s n , si and 9 n by 

t + u = 2M 2 + M 2 K + si-s n , 

v = -2a n Xcos6 n , (2.7) 



where 



a v = (l-4M 2 /s w )2 
X 



i 



^ 2 (M 2 K , S7T , Sl ) , 
A(x, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) . (2.8) 



In addition we define 



E = M 2 K + M 2 . (2.9) 



2.2 Matrix elements and decay rates 

The matrix element for K + — > 7r + 7r~/ + z/; is 



G 



T = ffiXPvh^ - IsWPiW* - A") , (2.10) 



where@ 



h = < 7T + ( Pl )n-(p 2 )ont | 1^(0) | K + (p) >; I = V,A , 
A = -^[P^ + Q^ + L^] , (2.11) 

M K 



1 In order to agree with the notation used by Pais and Treiman |l4| and by Rosselet et al. [|Tl| 
we have changed the previous convention M m in the definition of the anomaly form factor H. 
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and eoi23 = 1- The matrix elements for the other channels ( |2.2| , |2.3| ) may be obtained 
from (|2.10| , |2.11|) by isospin symmetry, see below. 

The form factors F, G, R and H are analytic functions of the variables s w , t and 
u. The partial decay rate for Q2.1|) is given by 



dY 



2M k (2tt) 8 



2 , T ,d 3 p 1 d 3 p 2 d 3 pi d 3 p v 

^ \T\ S (p-P-L)_ 



2p° 1 2p\ 2p Q t 2p° ' 



spins 



(2.12) 



The quantity 



spins 



T | 2 is a Lorentz invariant quadratic form in F, G, i? and 
H . All scalar products can be expressed in terms of the 5 independent variables 
Sn,si,9 n , 9i and 0, such that 



£ | T | 2 = 2G% | Ks | 2 M^ 2 J 5 K,^,^,^ 



(2.13) 



spins 



Carrying out the integrations over the remaining 4-3 — 5 = 7 variables in Q2.12 ) 
gives ||T3| 

dT 5 = G% | V us | 2 N(s n ,s l )J 5 (s n ,s h 9 7T ,9 h (j))ds 7T dsid(cos9 7T )d(cos9 l )d(j) , 
N(s„, Sl ) = (l-z l )a^X/(2 13 7 r 6 M 5 K ) , (2.14) 



where 



J5 = 2(1 — zi) h + h cos 29i + J 3 sin 2 Q\ ■ cos 20 + J 4 sin 29i ■ cos 1 
+ I 5 sin #j • cos 4> + Iq cos #z + I 7 sin • sin + J 8 sin 26 1 / • sin 
+ ig sin 2 61 ■ sin 2 



with 



h 

h 

h 

h 
h 
h 
h 

U 



i {(1 + ^)l^i| 2 + ~(3 + (|F 2 | 2 + |F 3 | 2 ) sin 2 6* + 2z,|F 4 | 2 } 
-j(l-^){l^ii 2 -5(l^i 2 + |F 3 | 2 )sin 2 ^} , 
-±(1 - z t ) {\F 2 \ 2 - \F 3 \ 2 } sin 2 9« , 

\(l-z{) Re(F*F 2 )sm9 n , 

- { Re(F*F 3 ) + m Re(F 4 *F 2 )} sin^ , 

- { Re(F 2 *F 3 ) sin 2 6 V - z x Re(F*F 4 )| , 

- { Im(F*F 2 ) + zi Im(F 4 *F 3 )} sin 9 n , 

Im(F*F 3 )sin^ , 

- (l-Zi) Im(F 2 *F 3 )sin 2 ^ , 



(2.15) 
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and 



Ft = X ■ F + a 7T {PL) cos 6 n ■ G 
F 2 = cr^ (s^si^^G , 

F 3 = a w X (s n si) lf ~ - — 



Ml 



F 4 = -{PL)F — siR — a n X cos 9 n ■ G . (2.16) 

The definition of Fx,...,F4 corresponds to the combinations used by Pais and 
Treiman |TJ| (the different sign in the terms ~ PL is due to our use of the metric 



diag(H )). The form factors I±, . . . , Jg agree with the expressions given in [14] 



2.3 Isospin decomposition 



The K14 decays (|2.2|) and (|2.3| ) involve the same form factors as displayed in Eq. 
( [2. 11|) . We denote by A 00 and A _ the current matrix elements of the processes 
( p.l|) -( p73|) . They are related by isospin symmetry!, 

A + - = ^=-A 00 . (2.17) 

This relation also holds for the individual form factors, which may be decomposed 
into a symmetric and an antisymmetric part under t <-> u (pi «-> pz)- Because of 
Bose symmetry and of the AI = | rule of the relevant weak currents, one has 

(F,G, R, H) 00 = —(F + ,G~,R + ,H ) + _ , 

(F, G, R, #)o- = V2(F~,G + ,R-,H + ) + _ , (2.18) 

where 

F±_ = -[F(s n ,t,u)±F{ Sir ,u,t)] , (2.19) 

and F(s n ,t,u) is defined in Eq. (|2.11|) . G ± ,R ± and are defined similarly. The 
isospin relation for the decay rates is 

T(K + -> ir+ir-l+Vi) = ^T(K° -> + 2T(K + -> 7r°7r°Z + ^) . (2.20) 

2.4 Partial wave expansion 

The form factors may be written in a partial wave expansion in the variable 8 n . We 



consider a definite isospin Tin state. Suppressing isospin indices, one has |TS 

-n_P_ 
X 



00 a PL 

Pi{cos 9 w )fi !L p- cos 9 W G 

1=0 



2 We use the Condon-Shortley phase conventions. Notice that we evaluate matrix elements and 
decay rates for K° - they differ from the corresponding A'i-quantities by a normalization factor. 
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G = £p/(cos^ 



1=1 

oo 



where 



1=0 A 

oo 

H = J2 P K cosd -) h i > ( 2 - 21 ) 
z=o 

P/(*) = -^P(z) . (2.22) 

The partial wave amplitudes fi,gi,ri and Zty depend on s n and s;. Their phase 
coincides with the phase shifts 8f in elastic 7T7t scattering (angular momentum Z, 
isospin J). More precisely, the quantities 

e~<X 2l , 

e-^+iX 2m ; Z = 0, 1, . . . ; X = f,g, r, h , (2.23) 

are real in the physical region of Ku decay (in our overall phase convention). The 
form factors Pi and P4 therefore have a simple expansion, 

1 

p 4 = -J2 p i( cosd *)( pL fi + s i r i)- ( 2 - 24 ) 



3 Theory 

The theoretical predictions of K14 form factors have a long history which started 
in the sixties with the current algebra evaluation of P, G, R and H. For an early 
review of the subject and for references to work prior to CHPT we refer the reader 
to |L6| (see also ||17|| ). Here we concentrate on the evaluation of the form factors in 



the framework of CHPT |, g 



3.1 Tree level 



The chiral representation of the form factors at leading order was originally given 
by Weinberg 0, 

P = G = = 3.74 , 

M K ( SjL + v \ 

H = . (3.1) 
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The next-to-leading order corrections are displayed below, and the later sections 
contain an estimate of yet higher order contributions. Here we note that the total 
decay rates which follow from Eq. (|3.1|) are typically a factor of two (ore more) 
below the data. As an example, consider the channel K + — ► 7r + 7r~ e + u e . Using 
11), the total decay rate becomes^ 1297 sec -1 , whereas the experimental value is 



3160±140 sec" 1 p] 



3.2 The form factors at one-loop 

In Ref. H m, the form factors F, G and H have been evaluated in CHPT at order 
p 4 (see also [19|). We complement these works with the evaluation of R at the same 
order. Below we display the result of our calculation, referring the reader to the 
above references and to available reviews || for the details of the methods used. 

In order to make this article reasonably self-contained, we display the result of 
all four form factors. The result for F may be written in the form 



F(s nj t, u) = jl + ±(U F + P F + C F ) + 0(E 4 ) 



(3.2) 



The contribution UF(s n , t, u) denotes the unitarity correction generated by the one- 
loop graphs which appear at order E A in the low-energy expansion. It has the form 



U F {s w ,t,u) = A (s 7r ) + Ap{t) + B(t, u) 



(3.3) 



with 

A (s 7r ) 
A F {t) 



3 s M 2 



16 
1 

+ r 



B(t,u) 



[UM 2 K + UMl - l9t)J r Kn (t) + (2M| + 2Ml - 3t) J'^t) 



(3M 2 K - 7M 2 + 5t)K Kn (t) + {Ml - 5M 2 + 3t)K nK {t) 



9(L K „(t) + L vK (t)) + (3M| - 3M 2 - M)(M r K „(t) + Mfc(t)) 
(M 2 K + M*-t)r K7r (t)-(t~u). 



(3.4) 



The loop integrals J^(s ir ), . . . which occur in these expressions are listed in appendix 
[A]. The functions J P q and M PQ depend on the scale /i at which the loops are 
renormalized. The scale drops out in the expression for the full amplitude (see 
below). 

The imaginary part of F 7 7 2 Ao(s 7r ) contains the 1 = 0, 5-wave nir phase shift 



S° (s n ) = (32tcF 2 )~ 1 (2s 7T - M 2 )a„ + 0(E 4 



(3.5) 



3 If not stated otherwise, we use = 93.2 MeV, \V US \ = 0.22 and (M^,M K ) = (139.6,493.6) 
MeV, (135,493.6) MeV and (137,497.7) Mev for the decays (|2j|) , (|J|) and (§J), respectively. 
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as well as contributions from KK and r]f] intermediate states. The functions Ap(t) 
and B(t,u) are real in the physical region. 

The contribution Pp(s 7r ,t, u) is a polynomial in s n ,t,u obtained from the tree 
graphs at order E A . We find 



P F (s n ,t,u) = ^2p itF (s n ,t,u)L r i , 



(3.6) 



i=i 



where 



Pl,F 
P2,F 
P3,F 
Pi,F 
P5,F 
P9,F 



32(s n -2M 2 ) , 
8(M 2 K + S7T - Sl ) , 
4(M| - 3M 2 + 2s w - t) 
32M 2 , 
AM 2 , 
2s, . 



(3.7) 



The remaining coefficients p^p are zero. The quantities L\ denote the renormalized 
coupling constants which parametrize the effective lagrangian at order E A ||. Their 
scale dependence is 

L r M=L r M + -^- 2 \n^ . (3.8) 

Observable quantities are independent of the scale /i, once the loop contributions 
are included. The coefficients Tj are displayed in table 1, together with the value || 
of the couplings L\ at \i = M p . 

The contributions Cp contain logarithmic terms, independent of s n ,t and u: 



C f 



256tt 2 



M 2 M 2 M 2 



LR ~ fi 2 



The corresponding decomposition of the form factor G is 



G(s n ,t,u) 
U G (s w ,t,u) 



M, 



K 



1 + -^(U G + P G + C G ) 



V2F, 

A 1 (s n )+A G (t) + B(t,u) 



0(E 4 



(3.9) 



(3.10) 



with 



= -[(2M 2 K + 2M 2 + 3t)J r K7T (t)-(2M 2 K + 2M 2 -3t)J r vK (t) 



+ 



{-3M 2 K + 7M 2 - 5t)K Kv {t) + (-M 2 K + 5M 2 - 3t)K r)K {t) 



L Kn {t) + L vK (t) - (M 2 K — M 2 + t)(M r Kw (t) + M r vK (t)) 



(3.11) 
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Table 1: Phenomenological values and source for the renormalized coupling con- 
stants Ll(Mp) according to Ref. ||. The quantities Tj in the fourth column deter- 
mine the scale dependence of the L\ (/i) according to Eq. ( |3.8| ). L r u and L r 12 are not 
directly accessible to experiment. 



i 


Ul(M p ) x 10 3 


source 


T t 


1 


0.7 ± 0.5 


7T7T D-waves, Zweij 


I rule 


3/32 


2 


1.3 ± 0.7 


mi .D-waves 




3/16 


3 


-4.4 ± 2.5 


7i7i D-waves, Zweij 


I rule 





4 


-0.3 ± 0.5 


Zweig rule 




1/8 


5 


1.4 ± 0.5 


Fr '■ F T 




3/8 


6 


-0.2 ± 0.3 


Zweig rule 




11/144 


7 


-0.4 ± 0.2 


Gell-Mann-Okubo 







8 


0.9 ± 0.3 


M K0 -M K+ ,L 5 , 




5/48 






(2m s -m u - m d ) : 


(m d - m u ) 




9 


6.9 ± 0.7 


< r 2 > 7r 

em 




1/4 


10 


-5.5 ±0.7 


71 — > eZ/7 




-1/4 


11 
12 






-1/8 
5/24 



The imaginary part of F n 2 Ai(s 7r ) contains the 1 = 1, P-wave phase shift 

SIM = {^txFI)- 1 s^I + 0(E A ) . (3.12) 

as well as contributions from KK intermediate states. The function Aq is real in 
the physical region. 
The polynomials 

9 

^ = EftA.*.«W ( 3 - 13 ) 
1=1 

are 

P2,G = S(t — U) , 

p 3 ,G = 4(t-M|-M 2 ) , 
p 5 ,G = 4M 2 , 
p 9<G = 2si , 

(3.14) 

The remaining p it c vanish. The logarithms contained in Cq are 

C G = -C F . (3.15) 

The form factor R contains a pole part Z(s n , t, u) / (si — M K ) and a regular piece 
Q. [Since the axial current acts as an interpolating field for a kaon, the residue of 
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the pole part is related to the KK — ■> ittt amplitude in the standard manner.] We 
write 



R 
I 



M 



K 



2V2F W [si-Mjc 



+ Q + 0(E 4 



B I + —(U I + P I + Ci), J = Z,Q 



(3.16) 



The Born terms Bj are 



So 



1. 



(3.17) 



For the loop corrections t/j, and Cj we find for the residue Z 



U z = s w A (s w ) + i/AxK) - -M 2 K M 2 J r m ( S7T ) 

+ 1 [ll( Sw - z/) 2 - 20S(^ - !/) + 12S 2 ] J^(t) 
+ l[ 3 (^-^)-2£] 2 J^(t) 



1 

I' 



, _( Sw + z/ )2j^( M ) 



+ -(M 2 K -M 2 )[5( S7T -u)-6^]K K7T (t) 
+ \{M 2 K -Ml)[^-v)-2T]K nK {t) 



2s^{v + 4S) - 3si + u 2 - 1QM 2 M 2 K M r Kn (t) + M r vK (t) 



and 



- -(3s 7r + u-2E)(L r , K (t) + L K „(t)) 



(3.18) 



with 



P z (s n ,t,u) 



^2p i ,z(s 7 :,t,u)L r i 



(3.19) 



Pl,Z 
P2,Z 
P3,Z 
Pl,Z 
P5,Z 

P6,Z 
P8,Z 



32(s n -2M 2 K )(s n -2M 2 ] 



(si + v 2 ) 



-2 

32 

4 



2(i/ + 4£K - 54 - i/ 2 - IQM 2 K M 2 



Es n - AM 2 K M 2 



S7T + v)Z- 8M 2 K M 2 
2 11 #2 



128M 2 K M n . 
UM 2 K Ml . 



(3.20) 
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The remaining p^z vanish. Finally, the logarithms in Cz are 

C z = 



Ml - Ml 



M 2 M 2 M 2 

3Ml In - 2M 2 K In ™£ - M 2 In ^ 



1287T 2 

The nonpole part Q receives the following one-loop contributions: 



and 



Ur 



M 2 K -si 
32 



{uj r K „(t) + 8J r K7r (u) + 3J r vK (t)} 



- ~{5{ S7r -v) + 5{M 2 K - Sl )-6E)K Kn {t) 

o 

l(3( S7T -u) + 3(M 2 K - Sl ) - 2Y)K nK {t) 

o 

9 

- -(L vK (t) + L Kv (t)) 

+ §(4(1/ + 2M, 2 ) - 3(M| - sOX^U*) + ^W), 



(3.21) 



(3.22) 



with 



P Q (s n ,t,u) 



$3 Pi,Q(sir, 



i=i 



32(^-2M 2 ) , 



p 2iQ = 8(M 2 K - Sl ) , 

p 3tQ = 2(4(s n -2Ml) + M 2 K - Sl ) 

p 4iQ = 32Ml , 

p 9iQ = 2[(s 7r + u)-(M] c -s l ) 

The remaining p^g vanish. Finally, the logarithms in Cq are 

1 



Cc 



1287T 2 



/If 2 M 2 M 21 

5 M 2 In ±% - 2M 2 K In ^ - 3M* In ^ 



(3.23) 



(3.24) 



(3.25) 



The first nonvanishing contribution in the chiral expansion of the form factor H 
is due to the chiral anomaly ||21|| . The prediction is @ 



H 



V2M 3 



K 



87T 2 Fi 



-2.66 



(3.26) 



in excellent agreement with the experimental value [TT[ H = —2.68 ± 0.68. The 
next-to-leading order corrections to H have also been calculated PD| . If the new 
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low-energy parameters are estimated using the vector-mesons only, these corrections 
are small. 

The results for F, G and R must satisfy two nontrivial constraints: i) Unitarity 
requires that F, G and R contain, in the physical region 4M^ < < [Mk — mi) 2 , 
imaginary parts governed by S- and P-wave hit scattering [these imaginary parts 
are contained in the functions A (s vr ), A 1 (s vr )]. ii) The scale dependence of the low- 
energy constants L\ must be compensated by the scale dependence of U f,g,z,q and 
Cf,g,z,q for all values of s n , t, u, M%, M\. [Since we work at order E A in the chiral 
expansion, the meson masses appearing in the above expressions satisfy the Gell- 
Mann-Okubo mass formula.] We have checked that these constraints are satisfied. 

It is one of the aims of this article to determine the low-energy constants L\, L r 2 
and L3 from experimental data on K + — > 7r + 7r~e + z/ e decays and on tttt threshold 
parameters. In Ref. 0, ||, the above one-loop expressions have been used for this 
purpose. Because the one-loop contributions are rather large, the values of the 
Lj's so extracted suffer from substantial uncertainties. In the following section, we 
therefore first estimate the effects from higher orders in the chiral expansion, using 
then this improved representation for the form factors in a comparison with the 
data. 

4 Beyond one-loop 

To investigate the importance of higher order terms, we employ the method devel- 
oped in Ref. |]12[ . It amounts to writing a dispersive representation of the partial 
wave amplitudes, fixing the subtraction constants using chiral perturbation theory. 
Here, we estimate the higher order terms in the S-wave projection of the amplitude 

F u ' 

= (47rX)- 1 | dfiF^Mi) , (4-1) 
because this form factor plays a decisive role in the determination of L"[, U 2 and L 3 , 



and it is influenced by S"-wave tttt scattering which is known [22, |2j| to produce 
substantial corrections. 

4.1 Analytic properties of f(s w ,si) 
Only the crossing even part 

F+ = XF + + a w (PL) cose n ■ G~ (4.2) 



contributes in the projection ( f4.1|) . The partial wave / has the following analytic 
properties: 

1. At fixed si, it is analytic in the complex s^-plane, cut along the real axis for 
Re s w > AMI and Re s w < 0. 

2. In the interval < s w < 4M^, it is real. 
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3. In 4M 2 < s n < 16MJ, its phase coincides with the isospin zero S-wave phase 
5q in elastic tttt scattering, 



f + = e 2 <f. ,f ± = f(s v ±ie, Sl ) 



(4.3) 



The proof of these properties is standard pi . Here we only note that the presence 
of the cut for s n < follows from the relations 



t = M* + M *" + Sl - S * 



— a w X cos 9^ 



t(cos6 n = -1, < 0) > (M K + M w ) 2 . 



(4.4) 



Since F + and G have cuts at t > (M K + M n ) 2 [see e.g. Eqs. |3~^ - |3.11|1 , the claim 
is proven. 

4.2 Unitarization 

We introduce the Omnes-function 



fl(s n ) = exp 



ds <J°(f 



71 S S 



(4.5) 



where A will be chosen of the order of 1 GeV below. According to (|4.3|) , multiplica- 
tion by f2 removes the cut in / for AM 2 < s n < 16M 2 . Consider now 

f = h + fn , (4-6) 
where /l(/r) has only the left-hand (right-hand) cut, and introduce 

v = JTV - f L ) . (4.7) 
Then t> has only a right-hand cut, and we may represent it in a dispersive way, 

si r°° dslmQ-^f - f L ) 



V = V o + ViS n + 



7T J AMI S 2 



s - S v 



(4- 



We expect the contributions from the integral beyond lGeV 2 to be small. Further- 
more, fi -1 / is approximately real between 16M 2 and lGeV 2 , as a result of which 
one has 

v = vq + vxSt, - — I . (4.9) 



71 JAM* S z S — S n 

For given vo, v i, /l and f2, the form factor / is finally obtained from 

f = h + n v . (4.io) 

The behaviour of Jl at s n — > is governed by the large |t|-behaviour of F + and 
G~, see ( |4.4|) . Therefore, instead of using CHPT to model /l, we approximate the 
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left-hand cut by resonance exchange. To pin down the subtraction constants vq and 
V\ , we require that the threshold expansion of / and / CHPT agree up to and including 
terms of order 0(E 2 ). For a specific choice of Jl, this fixes vq,vi in terms of the 
quantities which occur in the one- loop representation of F + and G~ . In the case 
where fi = 0, / has then a particularly simple form at si = 0, 



f(s n ,s l = 0)\f L=0 



v 



tt(v -\ 
Mk 

+2(M 2 



1.05 + 



F 2 



64M 2 L\ + 8M 2 K L r 2 



K 



8M 2 )L 3 + ^(M 2 K 



4M 2 )(4L r 2 + L 3 



M, 



K 



V2F n 



0.38 + - 



32L\ + 8L r 2 + 10L 3 



3 AM 2 V 23 



(4.11) 



We relegate the details of the calculation of Jl, Vq and of V\ to appendix |C[ 

In the partial wave /, the effects of the final-state interactions are substantial, 
because they are related to the 1 = 0, 5-wave tttt phase shift. On the other hand, 
for the leading partial wave in G + = ge l5p + • ■ •, these effects are reduced, because 
the phase S p is small at low energies. We find it more difficult to assess an estimate 
for the higher order corrections in this case - we come back to this point in the 
following section. 



5 Determination of L\,L<i and L3 

Here we determine the low-energy constants L\,U 2 and L 3 from experimental data 
on K + — > 7r + 7r~e + z/ e decays and on tttt — ► tttt threshold parameters, using the 
improved S'-wave amplitude / set up above. We are aware that our results will not 
be the last word: future kaon facilities like DA$NE [2(J will allow a more refined 
comparison of the chiral representation with the data. Nevertheless, we believe that 
it is instructive to see what one has to expect from higher order contributions. A 
comparison with earlier work |J will be provided at the end of this section. 

5.1 The data 

Experimentally the study of decays is dominated by the work of Rosselet et al. 
JllJ which measures K + — ► 7r + 7r~e + z/ e with good statistics. The total decay rate, 
the absolute value of the form factors F, G and of H and the difference of the phases 
Sq — 8{ were determined by use of 



F = 




+ f p e i5p cos 9 W + D — wave 


G = 


ge i5p 


+ D — wave , 


H = 


he iSh 


+ D — wave . 



(5.1) 
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The form factor f p was found to be compatible with zero and hence set equal to 
zero when the final value for g was derived. No s^, si dependence of the ratios gj f s 
and h/f s was seen. Parametrizing f s in the form 



2 



then gives 



9 
h 



/ s (0)(l + A /g 2 ) , 
(s^- AMI)/ AMI 



s(0)(l±A 9 g 2 
h(0)(l + X h q 2 



(5.2) 



(5.3) 



with \j = Aq = Xh- Rosselet et al. found |TT 



fs(0) 
9(0) 
h(0) 

A/ 



5.59 ±0.14 
4.77 ±0.27 
-2.68 ±0.68 
0.08 ±0.02 



(5.4) 



where we have used 
imental numbers 



\V US \ = 0.22 in transcribing their results. Notice that the exper- 



have been obtained in ]TT under assumptions which are in 
conflict with our theoretical formulae, like absence of higher waves, si independence 
of the form factors and equality of the slopes, \f = \ g = A^. It would of course be 
desirable to analyze forthcoming data without any additional assumptions. 



The total decay rate is flT 



T Ke4 = (3.26±0.15)10 3 sec 



-i 



(5.5) 



In H it has been observed that, using as input the central values (|5.4| ), one obtains 
^K e4 = 2.94- 10 3 sec -1 , which disagrees with ( |5.5[ ). On the other hand, the value (|5.5|) 
was used in Ref. 



Ill to normalize the form factors in 



We do not understand 



the origin of this contradiction. 

The 7T7T threshold data used below are taken from Ref. 
table 2, column 6. 



25]. We display them in 



5.2 The fits 

In the following, we perform various fits to / s (0), Xf,g(0) and to the 7T7T threshold 
parameters listed in table 2. We introduce for this purpose the quantities 



f(Sn,Sl) 

9(s n ,si) 



(47rX)- 1 | dnF 1 (s 7r ,t,s l 

3 r 

— I <ifisin 2 6 l 7r G(s 7r ,t, si) 

<57T J 



\f(Sn,Sl)\ 



(5.6) 



16 



Table 2: Results of fits with one- loop and unitarized form factors, respectively. The 
errors quoted for the L['s are statistical only. The L 4 r are given in units of at 
the scale \x = M p , the scattering lengths a{ and the slopes b\ in appropriate powers 
oiM n +. 





K e 4 data alone 


K e4: and 


Tin data 


experiment 




one-loop 


unitarized 


one-loop 


unitarized 


13 


L\ 


0.65 ±0.27 


0.36 ±0.26 


0.60 ±0.24 


0.37 ±0.23 




L r 2 


1.63 ±0.28 


1.35 ±0.27 


1.50 ±0.23 


1.35 ±0.23 




L 3 


-3.4 ± 1.0 


-3.4 ±1.0 


-3.3 ±0.86 


-3.5 ±0.85 




a° 


0.20 


0.20 


0.20 


0.20 


0.26 ±0.05 


b°o 


0.26 


0.25 


0.26 


0.25 


0.25 ±0.03 


-10 a 2 


0.40 


0.41 


0.40 


0.41 


0.28 ±0.12 


-10 bl 


0.67 


0.72 


0.68 


0.72 


0.82 ±0.08 


10a\ 


0.36 


0.37 


0.36 


0.37 


0.38 ±0.02 


10 2 b\ 


0.44 


0.47 


0.43 


0.48 




10 2 a° 


0.22 


0.18 


0.21 


0.18 


0.17 ±0.03 


10 3 a 2 


0.39 


0.21 


0.37 


0.20 


0.13 ±0.3 


X 2 /N DOF 


0/0 


0/0 


8.8/7 


4.9/7 





where the factor 3/2 sin 2 9 appears because G is expanded in derivatives of Legendre 
polynomials. Below, we confront [/ s (0), g(0)] with [/(4M 2 ,S;)], g(4M%, si)], which 
depend on sj. Furthermore, we compare the slope Xf with 

which depends on both s n and sj. Below we use these dependences to estimate 
systematic uncertainties in the determination of the low-energy couplings. [In future 
high statistics experiments, the ^-dependence of the form factors will presumably 
be resolved. It will be easy to adapt the procedure to this case.] 

We have used MINUIT |26] to perform the fits. The results for the choice si = 
0, s n = 4.4M 2 are given in table 2. In the columns denoted by "one-loop", we have 
evaluated f,g and A/ from the one- loop representation given aboveQ. In the fit 
with the unitarized form factor (columns 3 and 5), we have evaluated / from Eqs. 
QOOyCH] ), inserting in the Omnes function the parametrization of the txtx S-wave 



phase shift proposed by Schenk [27|, solution B]. For the form factor G, we have 
again used the one-loop representation. The statistical errors quoted for the L^s 
are the ones generated by the procedure MINOS in MINUIT and correspond to an 
increase of x 2 by one unit. 

A few remarks are in order at this place. 



We always use for L\,...,L\ the values quoted in table 1. 
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1. 



It is seen that the overall description of the tttt scattering data is better using 
the unitarized form factors, in particular so for the D-wave scattering lengths. 



The errors quoted do not give account of the fact that the simultaneous deter- 
mination of the three constants produces a strong correlation between them. 
To illustrate this point we note that, while the values of the Lj's in column 4 
and 5 are apparently consistent with each other within one error bar, the x 2 
in column 5 increases from 4.9 to 30.7 if the Lj's from column 4 are used in 
the evaluation of x 2 i n column 5. (For a discussion about the interpretation 
of the errors see 



The low-energy constants li, l 2 which occur in SU (2) L x SU{2) R analyses may 
be evaluated from a given set of L\ , L r 2 and L3 [|3] . Their value changes in 
a significant way by using the unitarized amplitude instead of the one-loop 
formulae: the values for (Zi, I2) in column 4 and 5 are (—0.5 ± 0.88, 6.4 ± 0.44) 
and (-1.7 ± 0.85, 6.0 ± 0.4), respectively. 



4. L\, U 2 and L 3 are related to tttt phase-shifts through sum rules [^8|, |29|| . In 
principle, one should take these constraints into account as well!. We do not 
consider them here, because we find it very difficult to assess a reliable error 
for the integrals over the total 7T7T cross sections which occur in those relations. 

The statistical error in the data is only one source of the uncertainty in the low- 
energy constants, which are in addition affected by the ambiguities in the estimate 
of the higher order corrections. These systematic uncertainties have several sources: 

i) Higher order corrections to g have not been taken into account. 

ii) The determination of the contribution from the left-hand cut is not unique. 

iii) The quantities / and g depend on si, and A/ is a function of both si and s w . 

iv) The Omnes function depends on the elastic 7T7T phase shift and on the cutoff 
A used. 

We have considered carefully these effects. As for the first point, we have evalu- 
ated the higher orders in g in two ways: 



We define the quantity 



A - (4(0) ~^ (2) ) 2 
9 = P ' 



where g® is the CHPT prediction at leading-order. We then add Ag in 
quadrature to the experimental error in g(0) and redo the fit. This generates 
slightly larger errors than before. To illustrate, the entries (0.23, 0.23, 0.85) in 
column 5 in table 2 become (0.29, 0.28, 1.1). 



3 We thank B. Moussallam for pointing this out to us. 
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Table 3: Fits made with different theoretical input to determine the systematic 
uncertainties in L"[, L r 2 and L 3 . AL\ shows the difference to the value of L\ displayed 
in the second row. L\ are given in units of 10~ 3 at the scale \x = M p . The fits include 
K ei and txtx data. 





L\ ALi 


U 2 AL 2 


L 3 AL 3 


5 th column table 2 


0.37 


1.35 


-3.5 


full propagators 

h = o 

s w = 5.6M%,si = 
s n = AAMl Sl = M 2 K / 10 
txtx phase from 27, solution A] 
txtx phase from []27|, solution C] 


0.34 -0.03 
0.42 0.05 
0.38 0.01 
0.33 -0.04 
0.33 -0.04 
0.46 0.09 


1.08 -0.27 
1.30 -0.05 
1.35 0.00 
1.26 -0.09 

1.35 0.00 

1.36 0.01 


-2.8 -0.7 
-3.5 0.0 
-3.5 0.0 
-3.3 -0.2 
-3.5 0.0 
-3.5 0.0 



The main contribution to the one- loop correction in g stems from L 3 . On the 
other hand, the low-energy constants L\ are saturated by resonance exchanges 
whose contribution is evaluated in the limit of large resonance mass |3(| . We 



find it therefore resonable to estimate higher orders by using the complete 
resonance propagators. In order to evaluate the uncertainty induced by this 
correction, we made the fit including the full propagators in g and in /. This 
changes the central values for the Li's, see table 3. 

Concerning the effect of the left-hand cut, we estimate its uncertainty by dropping 
this piece altogether. The results of the fit in this case are again given in table 
3. It is seen that the L['s depend rather weakly on the presence of the left-hand 
cut. Figures 2 and 3 illustrate the effect of Jl in more detail. In Fig. 2 is shown 
the form factor / with and without the left-hand cut (dashed and dash-dotted line, 
respectively). In both cases, L[,L 2 , L 3 from ( |5.9| ) have been used. At large values of 
s n , the difference between the two evaluations of / is not negligible for the following 
reason. At threshold, the value of the form factor and of its slope cannot change 
much by construction-these quantities are constrained in our procedure by the low- 
energy expansion at one-loop order. On the other hand, the quadratic piece is 
unconstrained and receives contributions from fi through the dispersive integral. 
The corresponding change in / is less than 10% at = 10M^. In Fig. 3, we 
display the differential decay rate dT/ds n , obtained by using the form factor with 
and without /x,. The figure illustrates that the contribution from the quadratic 
piece in to observables is strongly suppressed by phase space and will therefore 
be difficult to observe. 

The dependence of the fits on si, s n and on the txtx phase shift used in the Omnes 
function Q is illustrated in the last four rows in table 3. Furthermore, changing 
the cutoff A = lGeV in Q to A = 0.8GeV induces small changes in the Lj's only. 
We conclude that a global fit to all the available data is rather stable against the 
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systematic uncertainties considered here. 

To finally give the best determinations of L\,U 2 and L 3 , we take the central 
values from the global fit displayed in table 2, column 5. For the corresponding 
errors, we take the ones generated by using the theoretical error bars for the higher 
orders in g, and find in this manner 

10 3 L^(M p ) = 0.4 ±0.3 , 

10 3 L r 2 {M p ) = 1.35 ±0.3 , (5.9) 
10 3 L 3 (M p ) = -3.5 ±1.1 . 

For SU(2)l x SU(2)r analyses it is useful to know the corresponding values for the 
constants l\ and I2, 

h = -1.7 ±1.0 , 

l 2 = 6.1 ±0.5 . ^- iUj 

The value and uncertainties in these couplings play a decisive role in a planned 
experiment JJlj] to measure the lifetime of tt + 7t~ atoms, which will provide a com- 
pletely independent measurement of the tttt scattering lengths |a° — cLq\- 

One motivation for the analysis in f|, [| was to test the large Nc prediction 
U 2 = 2L\. The above result shows that a small non-zero value is preferred. To obtain 
a clean error analysis, we have repeated the fitting procedure using the variables 



x l = 


jr 
L 2 


2L\-U , 


xi = 


jr 

^2 ) 




X3 = 


(LI- 


-2L[)/L 3 . 


and mi 




includine the theoretical error in G 



discussed above, and found 

Xi = (4.8 ±0.8) • 10~ 3 , 

X 3 = -0.17t8;£ . (5.11) 
The result is that the large A^^ prediction works remarkably well. 

5.3 Comparison with earlier work 

It is of interest to compare the present procedure to determine the low-energy con- 
stants Lj_, U 2 and L3 with the method used in 0. There are two main differences: 

1 . The definition of the slope A / and of the threshold value of the form factors / s , g 
chosen in differs from the one used here. These quantities have of course a 
unique meaning in principle - on the other hand, one may wish to approximate 
a particular experimental situation. The procedure used in was adapted to 
Ref . |IT| , whereas a slight variation of the method proposed here may be useful 



once the s/-dependence of the form factors has experimentally been resolved. 
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2. Higher order corrections are estimated in || in a rather crude manner. In the 
present approach, the final-state interactions in the 1 = 0, S-w&ve amplitude 
are instead taken into account, and higher order terms in g are estimated with 
resonance exchange. 

The main effect of these differences can be described as follows. The different 
slope and form factors used in || lead to slightly different central values for L\, L 2 
and L3 at one-loop order, whereas the errors turn out to be very similar in both 
cases. The higher order estimates in |J lead to the same central values with larger 
error bars, whereas the unitarization performed in the present work leads to different 
central values with slightly smaller error bars than before, see columns 2/3 and 4/5 in 
table 2. This effect can be easily understood by considering the simplified expression 
( f4.11| ), which shows how the Omnes function affects the influence of the Lj's and 



hence their value in the fit. 



5.4 Improvements 



As we mentioned at the beginning of this section, there is room for improvement in 
the above treatment, both on the theoretical and on the experimental side. Con- 
cerning the latter, one should determine in future experiments the form factors f s 



and g without additional assumptions |TT[ which are in contradiction with the chiral 



representation. It remains to be seen whether this can be achieved by comparing 
the data directly with a modified chiral representation. In the latter, the full S- 
and P- wave parts of Fi and F 2 could be inserted, using the chiral representation 
solely to describe the small background effects due to higher partial waves I > 2. 
To be more precise, one would take for R and H the one-loop chiral representation, 
whereas for G one writes 

G = g(s w ,s l )e iS '+AG , 
3 f 

AG = G CH PT / dfisin^GcHPT , (5.12) 



8tt 

and similarly for F. The unknown amplitudes g(s n , si), f s (s n , si) and the phases 
S p , S s would then be determined from the data. We have checked that, if the errors 
in the form factors determined in this manner can be reduced by e.g. a factor 3 with 
respect to the ones shown in (|5.4j) , one could pin down particular combinations of 
L"[, L 2 and L3 to considerably better precision than was shown above. This is true 
independently of an eventual improvement in the theoretical determination of the 
higher order corrections in the form factor G - which is a theoretical challenge in 
any case. 

6 Predictions 

Having determined the constants L\, L r 2 and L 3 , there are several predictions which 
we can make. Whereas the slope \ g was assumed to coincide with the slope Xj in 
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the final analysis of the data in Ref. |TTJ, these two quantities may differ in the 
chiral representation. Furthermore, our amplitudes allow us to evaluate partial and 
total decay rates. In this section, we consider the slope X g and the total rates. 



6.1 The slope X g 

We consider the form factor g introduced in ( |5.6| ) and determine its slope X g 

g{s^ si) = g(4M 2 n , sl)(l + \ g ( Sl )q 2 + 0(g 4 )) (6.1) 

from the one-loop expression for G. The result is A g (0) = 0.08. As the slope is 
a one- loop effect, higher order corrections may affect its value substantially. For 
this reason, we have evaluated X g also from the modified form factor obtained by 
using the complete resonance propagators (and the corresponding Li a), compare 
the discussion above. The change is A\ g = 0.025. We believe this to be a generous 
error estimate and obtain in this manner 

X g (0) = 0.08 ±0.025 . (6.2) 

The central value indeed agrees with the slope Xf in (|5.4j ). 

6.2 Total rates 

Once the leading partial waves / and g are known from e.g. K + — > 7r + 7r~e + z/ e decays, 
the chiral representation allows one to predict the remaining rates within rather 
small uncertainties. We illustrate the procedure for K + — > 7r°7r°e + z/ e . According 
to Eq. (|2.18|) , the relevant amplitude is determined by F + ,G~,R + and H~. The 
contribution from H is kinematically strongly suppressed and completely negligible 
in all total rates, whereas the contribution from R is negligible in the electron modes. 
Using the chiral representation of the amplitudes F + and G~, we find that the rate 
is reproduced to about 1%, if one neglects G~ altogether and uses only the leading 
partial wave in the remaining amplitude, ~ —Xf . From the measured |TT]] form 



factor / = 5.59(1 + 0.08g 2 ) we then find r^+_ )7r o 7r o e + l/e = 1625sec 1 . Finally, we 
estimate the error from 

Ar = {[T(f s (0)+^fs,X f )-T(f s (0),Xf)] 2 + 

[r(/ s (o), Xf + AXf) - r(/ s (o), A/)] 2 } 1/2 = oosec" 1 , 

(6.3) 

where Af s = 0.14, AXf = 0.02. The final result for the rate is shown in the row 
"final prediction" in table 5, where we have also listed the tree and the one-loop 
result, together with the experimental data. The evaluation of the remaining rates 
is done in a similar manner - see table 4 for the simplifications used and table 5 for 
the corresponding predictions. 
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Table 4: Approximations used to evaluate the total rates in table 5. Use of / 
/ CH p T ,3 = (7chpt reproduces the one-loop results in table 5 to about 1%. 



a) K + decays 





7r°7r°e + z/ e 


1X + 1X + 


7T 7T / U + Z/ M 




-xf 


Xf + a w (PL) cos 6 n g 


-Xf 


F 2 










F 3 











F 4 


{PL)f 


- {(PL)f + siRchpt + cos 9 w g} 


{{PL)f + Sl R CiiPT } 



b) K° decays. Shown are the amplitudes divided by v2. 





7T U 7T e + U e 


7T U 7T ll + U^ 


F 1 


XF~ HPT + a n (PL) cos O^g 


XF~ HPT + av(PL) cos O^g 


F 2 




a^s^SiY^g 


F 3 








F 4 


- {(PL)F~ npT + a v X cos 6^} 


- {(PL)F c - pT + siR CUPT + a^XcosO^g} 



We have assessed an uncertainty due to contributions from F~ HPT , R CHPT in the 
following manner, i) We have checked that the results barely change by using the 
tree level expression for -R C hpt instead of its one-loop representation. We conclude 
from this that the uncertainties in -R C hpt do not matter, ii) The uncertainty from 
F~ UPT is taken into account by adding to Ar in quadrature the change obtained by 
evaluating F~ HPT at L 3 = —3.5 + 1.1 = —2.4. iii) In K° decays, we have also added 
in quadrature the difference generated by evaluating the rate with M n = 135 MeV. 

The decay K° — > 7r°7r~e + z/ e has recently been measured |33| with considerably 
higher statistics than before f3"2fl . We display the result for the rate in the first 
column of table 5b. The quoted errors correspond to the errors in the branching 



ratio and do not include the uncertainty in the total decay rate quoted by the 
PDG [32|. Notice that the value for L 3 determined in 
-1 " 



should be multiplied with 



7 Summary and conclusion 

1. The matrix elements for Ku decays depend on four form factors F, G, H and 
R. This article contains the full expressions for these at order E 4 in CHPT, 
thus completing already published calculations [?|, H, § of F, G, H at this order. 
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Table 5: Total decay rates in sec -1 . To evaluate the rates at one-loop accuracy, 
we have used L\,U 2 and L3 from ( |5.9| ). The remaining low-energy constants are 
from table 1. The final predictions are evaluated with the amplitudes shown in 
table 4, using / = 5.59(1 + 0.08g 2 ),# = 4.77(1 + 0.08g 2 ). For the evaluation of the 
uncertainties in the rates see text. 



a) K + decays 





7T + 7T e + V e 


7r°7r°e + z/ e 


7T + 7r _ / U + ^ 1 




tree 


1297 


683 


155 


102 


one-loop 


2447 


1301 


288 


189 


final 


input 


1625 


333 


225 


prediction 




±90 


±15 


±11 


experiment 


3160 


1700 


1130 




@ 


±140 


±320 


±730 





b) K° decays 





7T°7T e + V e 




tree 


561 


55 


one-loop 


953 


94 


final 


917 


88 


prediction 


±170 


±22 


experiment 


998 




m 


±39 ± 43 





We have estimated higher order terms in the S"-wave amplitude of the form 
factor F by use of a dispersive representation, determining the subtraction 
constants in the standard manner |]12| from CHPT. This procedure puts earlier 



attempts to estimate these corrections on a more firm basis. 

Using the improved S- wave amplitude, we have determined L[, L r 2 and L 3 from 
K + — > 7r + 7r~e + z/ e decays and tctt threshold data. Unitarizing the amplitude 
affects the related SU (2) x SU (2) constant Zi in a significant manner. As 
a result of this, the D-wave scattering lengths are in better agreement with 



the values given by Petersen |25fl than was the case before 0. All in all, a 



remarkably good agreement with K e \ and txix data is obtained. 

K14 decays may be used to test the large- Nc prediction L r 2 = 2L\ f|, |[ . Using 
the improved representation of the amplitudes, we have confirmed the earlier 
finding: The large- Nq rule works at the one standard deviation level for 
this combination of the constants. 
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5. The above determination of L\, U 2 an d -^3 w * n presumably be even more reli- 



able, once high statistics data from kaon facilities like DA$NE [20] will become 
available. 

6. We also predict the slope \ g of the form factor G and total decay rates, see 
Eq. (O) and table 5. 



We have made some effort to find out whether any of the decays could 
serve to determine some of the other low-energy constants which occur in the 
amplitude. We believe that it will be very difficult to pin down any of these 
(in particular L£) to better precision than already known, because the higher 
order corrections tend to wash out their effect. 

Finally, we would like to recall that the determination of the low-energy con- 
stants from K14 decays or the prediction of the total rates is not the only issue: 
these decays are in addition the only known source for a precise determina- 
tion of the isoscalar 5-wave phase shift near threshold. The possibilities 
to determine those in future high statistics experiments are presently under 
investigation |35| . 
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A Loop integrals 

In this appendix we define the loop integrals used in the text. We consider a loop 
with two masses, M and m. All needed functions can be given in terms of the 
subtracted scalar integral J(t) = J(t) — J(0) evaluated in four dimensions, 

. r d d p 1 
W = ~ % 1 (2n) d ((p + k) 2 - M 2 )(p 2 - m 2 ) ' ( ' ' 

with t = k 2 . The functions used in the text are then : 
-, . 1 f 1 , , M 2 -tx(l-x) -Ax 



16tt 2 Jo M 2 -Ax 

If A m 2 E, m 2 v/A, (t + V^) 2 - A 2 
2 H log — — ■ — — log —— log 



32tt 2 [ t M 2 A M 2 t to {t - v^) 2 - A 
r(t) = J(t)-2k, 

M r (t) = —{t-2YAJ(t) + ^J(t) + —^—-- 
w 12t 1 1 w 3t 2 W 288vr 2 6 

1 f_ M 2 m 2 , m 2 1 
E + 2—— log ■ 



96?r 2 t [ A M 2 J ' 

A 2 - 

= 

A = M 2 — m 2 , 
E = M 2 + m 2 , 

A = A(t,M 2 ,m 2 ) = (t + A) 2 -4tM 2 . (A.2) 

In the text these are used with subscripts, 

J ij {t) = J(t) with M = M i ,m = M j , (A.3) 

and similarly for the other symbols. The subtraction point dependent part is con- 
tained in the constant k 

1 M 2 log(f)-^log(^) 
32tt 2 M 2 -m 2 ' 1 ' 

where \i is the subtraction scale. 



B Resonance contribution to the form factors 



Below we display the contributions to the form factors F and G from resonance 
exchange (spin less than or equal to one, see also |36|, |37], |30[|). We quote them for 
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K + — > 7t + 7t~1 + i>i. The others can be derived using isospin relations ( |2.18| ). These 
contributions have been used both to provide a reasonable approximation of the 
left-hand cut, and to estimate higher order corrections in g. 



VECTORS 

• t-channel 



Fy 



M K G, 



x 



s^-channel 



2y/2F% M 2 - 1 
~F v (t -u+ 2si) + G v {t - u - 3s„ - Sl + M 2 K + 8M*) 
M K G V 1 
2x/2F3 M 2 - t X 

F V {M 2 K + si- SlT ) + G v (t -u + M 2 K + s n -si) 
M K G V 1 



V 



[F v - 2G V ] (t - u) 



2y/2F* M v - s, 

M K Gy 1 

2^2Fl M v -s w X 

(F v - 2G v )( Sl - M\) + (F v + 2GVK 



(B.1) 



(B-2) 



SCALARS 
• t-channel 



s^-channel 



\/2M 



K 



F* Ml-t 



c\(Ml + Ml -t)- c d c m (M 2 K + Ml) 



2y/2M- 1 



K 



Fi 



Ml - s w 



x 



c 2 d {2Ml - s n ) - 2c d c m Ml 



G 



5 







The values used for the couplings Gy, F v and c d , c m are |30| 

F v = 154 MeV , 
G v = 69 MeV , 
c d = 32 MeV , 
c m = 42 MeV , 



(B-3) 



(B.4) 



(B.5) 
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while for the masses we used 



M v = 770 MeV , 

M s = 985 MeV. (B.6) 

C Evaluation of /l,^o an d v\ 

fi is calculated from 

h = ^ / dQ + F+) + ^jp cos^(G y + . (C.l) 

Only the t-channel contributes to The s-channel has only singularities on the 
right-hand cut. The quantities in Eq. ( |C.1|) are defined as: 

= ^(45 + Ji v 1 s(^«)) , 
Gy, s = \{Gv,s-G v , s (t~u)) , (C.2) 

analogously to ( |2.19| ), see appendix B for Gy S ,Fy S . To evaluate v and V\, we 
impose that the unitarized amplitude f = fi + VLv matches the chiral one-loop 
representation / CH p T at the threshold s n = 4M^. We write 



M. 



/cHPT — (/cHPT + /cHPT + 0(E )) , 

/cHPT = 1 5 (C-3) 

with obvious notation, and have 

/&* + = + (1 + A) K } + ^ + «iV + 0(i? 4 )) , (C.4) 

where 

r A2 ds 5 °(s) 



7T J4M5 S S 



Is -Ml 



S °o = -g^f V 1 - 4 ^ , (C5) 
and where 5g is a quantity of order E 2 . The quantity 1 + A is the expansion of Q 

(k) 

in CHPT to the required order. The v\ are obtained by equating the threshold 
expansion of the left-and right-hand side in ( |C.4|) . 
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Figure captions 



Fig. 1 Kinematic variables for Ku decays. The angle 9 n is defined in T, 27T} 0i in 
and 4> in 

Fig. 2 The partial wave amplitude f(s w , S[ = 0). The dashed line shows /, evaluated 
according to Eqs. ( [4.10| , |C.1| ), with L\,U 2 and L 3 from (|5.9|). The dash-dotted 
line is evaluated with /l = according to (}4.11|) , using the same Lj's, whereas 
the solid line displays f s = 5.59(1 + 0.08g 2 ). 

Fig. 3 Differential rate dT/ds^ for K + — > 7r + 7r~e + z/ e decays in arbitrary units. The 
evaluation is done with Ff = F 2 = F 3 = F 4 ~ = 0, and F* = X f(s n , si = 
0), F 4 + = —PL/XF^. The input for the dashed (dash-dotted) line is the same 
as for the dashed (dash-dotted) line in Fig. 2. 
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